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■ ABSTRACT 

Context. The high-redshift (z = 3.1) blazar PKS 0537-286, belonging to the flat spectrum radio quasar blazar subclass, is one of the 
most luminous active galactic nuclei(AGN) in the Universe. Blazars are very suitable candidates for multiwavelength observations. 
Indeed, the relativistic beaming effect at work within the jet enhances their luminosity. This in turn allows the properties of the 
extragalactic jets, the powering central engine, and the surrounding environment to derived. 

Aims. Our aim is to present the results of a multifrequency campaign from the near-IR to hard X-ray energies on PKS 0537-286 and 
give insight into the physical environment where the radiation processes take place. 



o 

Methods. We observed the source at different epochs from 2006 to 2008 with INTEGRAL and Swift, and nearly simultaneously with 
ground-based optical telescopes. We also analyzed two archival spectra taken with XMM-Newton in 1999 and 2005. A comparative 
analysis of the results is performed. 

Results. The X-ray continuum of the blazar, as sampled by XMM, is described by a power law of index T = 1.2, modified by 
variable absorption at the soft X-rays, as found in other high-redshift QSOs. Modest X-ray continuum variability is found in the Swift 
^/-^ ' observations. The combined SwiftfBAT and SwiftfXKI spectrum is very hard (T = 1.3). This, together with the non simultaneous 

\Q | EGRET detection and the more recent non detection by Fermi-LAT, constrains the peak of the high-energy component robustly. 

, The optical/UV data, heavily affected by intervening Ly a absorption, indicate the presence of a bright thermal accretion disk that 

■ decreased in luminosity between 2006 and 2008. We infer from this a reduction of the BLR radius. When taking this into account, 

the 2006 and 2008 SEDs are compatible with a model based on synchrotron radiation and external inverse Compton scattering where 
the accretion-disk luminosity decreases between the 2 epochs by a factor 2, while the bulk Lorentz factor remains unchanged and the 
magnetic field changed only marginally. 
G\ . 

Key words, quasars: individual (PKS 0537-286) - galaxies: active - X-rays: galaxies - galaxies: high-redshift - radiation mecha- 
nisms: nonthermal 



^ ■ 1. Introduction sponsible for the low-energy hump generates the high-energy 

' hump through the inverse-Compton mechanism dSikora et alj 

Simultaneous multiwavelength studies of blazars offer a tool for [T^. | Dermer & Schlickeiseri[l993t iGhisellini & Madaul l l996). 

determining the properties of extragalactic jets, their powering Another mode i pre dicts that t he seed photo ns fo r the inverse- 

central engine, and the surrounding environment. The alignment Compton are external to the jet (ISikora et al.lll 994l:lDermer et all 

ofthe let with respect to theobserver's line of sight CAntonuccjl [199!. Alte rnative explanations are proton-initiated cascades 

[19931 l Urry & Padovani| [1995D makes blazars very luminous (M annheim 1993) or m-oton-svnchrotron emission from ultra- 



from radio frequencies to gamma-ray energies (lUlrich et al 



1997). Indeed, the beaming in these relativistic outflows (Rees 



relativistic protons (Mtic ke et alJl2003tlAharo nian 2000). 



19661) enhances the luminosity. While the jet propagates, part of High-redshift blazars are important tracers of the evo- 



its power is dissipated via cooling and acceleration processes of lution of super massive bla c k hol es (SMBH) in the early 
protons and leptons. The emitting particles, together with the rel- Universe. iBottcher & Dermerl d2002l) link the two blazar sub- 
ativistic beaming, give rise to the typical two broad hump struc- classes, namely flat spectrum radio quasars (FSRQ) and BL 
ture of the spectral energy distribution (SED). The low-energy the Lacertae (BL Lac) objects, which are the beamed counter- 
emission is understood and is believed to be produced through parts of low- and high-luminosity radio galaxies, respectively 
synchrotron radiation of relativistic electrons. In contrast, the JWall & Jackal 1 997b IWillottet alj|200ll) , through the gradual 
nature of the high-energy hump is still being debated. A possi- depletion of the circum-nuclear environment of SMBHs. In a 
ble explanation is that the same electron population that is re- process that is not understood well, the interaction of the ac- 

cretion disk and the SMBH drives the formation of jets. The 

Send offprint requests to: E. Bottacini, e-mail: eub@mpe.mpg . de blazar spectra therefore carry the information about the accre- 



2 



E. Bottacini et al.: PKS 0537-286, a multifrequency campaign 



tion and the circumnuclear environment, hence on the evolution 
of the environment of SMBHs over cosmic time. Indeed, a gen- 
eral characteristic of radio-loud quasars at high re dshift is an 
apparent absorption in excess of the Galactic value dFiore et al.l 
19981: [Reeves & Turner 2000). The nature of this X-ray spectral 
feature is still controversial. Obtaining precise spectra of blazars 
allows us to understand the physics at work in jets and derive 
information on the environment where the jet propagates. This 
in turn can confirm or deny evoluti onary scenarios of b lazars. 

PKS 0537-286, at z = 3.104 JWrightetalJll978h . is one 
of the most luminous kno wn high-redshift q uasars. First de- 
tected at radio frequencies dBolton et al.lll975 | ). it was observed 
at X-rays by the Einstein obs ervatory dZamorani et all 19811) 
and l ater studie d with ASCA dCappi et alj |1997t ISiebert et all 
1996), R OSAT (iFiore et al.1 fl998l) . XMM dReeves et alj 12001 
and Swift (Sambrun a et alj2007l) . These observations showed an 
extremely luminous quasar (L A = 10 47 erg s _1 in the 0.1-2 keV 
range) with a particularly hard spectrum (F = 1 measured by 
Swift/BAT) and a weak iron K emission l ine an d reflection fea- 
tures. Moreover. ISowards-Emmerd et al.l (|2004) identified PKS 
0537-286 as probable counterpar t of the EGRET source 3EG 
J053 1-2940 (lHartman et alJI 19991) . 

In this paper we report on a multiwavelength campaign on 
PKS 0537-286 ranging from the near-IR/optical to hard X-ray 
energies. In S. 2 we describe the observations of the 2006 and 
2008 multiwavelength campaigns; and in S. 3 the archival 2005 
XMM observations; in S. 4 we report the results of our multi- 
wavelenght observations and of the XMM archival data; in S. 4 
we discuss the results. 



Table 1. Log of UVOT observations. Flux density in units of 

10~ 14 erg cirT 2 s -1 A -1 



UTC start date 


exposure (ks) 


U-band flux density 


2006-10-27 05:55:49 


1.007 


(4.6 ± 1.0) 


2006-10-30 06:15:23 


1.090 


(5.9 ± 1.0) 


2006-10-31 06:21:27 


1.084 


(4.2 ± 1.0) 


2008-02-10 16:22:46 


0.484 


(5.6 ± 1.2) 


2008-02-12 16:24:40 


0.925 


(3.9 ± 1.1) 



Table 2. Log of Palomar observations. " statistical errors only. 



MJD 


UTC 


band 


Mag 


Err" 


Seein 


54033.40301 


2006-10-25 09:40:20 


B 


19.58 


0.02 


2"4 


54033.41492 


2006-10-25 09:57:29 


V 


18.90 


0.01 


2"0 


54034.42265 


2006-10-26 10:08:36 


B 


19.31 


0.03 


3"8 


54034.43479 


2006-10-26 10:26:05 


V 


18.86 


0.02 


2"5 


54035.41863 


2006-10-27 10:02:49 


B 


19.59 


0.03 


3"7 


54035.43105 


2006-10-27 10:20:42 


V 


18.99 


0.02 


3"1 


54036.41601 


2006-10-28 09:59:03 


B 


19.68 


0.02 


2"8 


54036.42706 


2006-10-28 10:14:57 


V 


18.86 


0.02 


2"5 


54038.40941 


2006-10-30 09:49:33 


B 


19.55 


0.03 


3"7 


54038.42128 


2006-10-30 10:06:38 


V 


18.96 


0.01 


3"4 


54039.40719 


2006-10-31 09:46:21 


B 


19.64 


0.02 


2"7 


54039.41966 


2006-10-31 10:04:18 


V 


18.89 


0.01 


2"3 


54040.40498 


2006-11-01 09:43:10 


B 


19.57 


0.02 


2"0 


54040.41697 


2006-11-01 10:00:26 


V 


18.67 


0.06 


1"7 



2. Observations 

2.1. The 2006 to 2008 multiwavelength campaign 

The source was monitored by INTEGRAL. Since the SwiftfBAT 
light curve over 9 months of exposure (lAiello et al.l [2008 a) 
showed a constant flux level, the INTEGRAL pointings were 
proposed in non-contiguous time intervals. ISGRI detected the 
source at 5.6 cr. In order to cover simultaneously the X-ray 
and UV ranges we requested Swift ToOs. These observations 
were scheduled on October 27"', 30 f \ 31 sf 2006 and further 
Swift/XRT observations were executed on February 10' /! , 12''' 
2008. Simultaneously to the Swift pointings REM, GROND, and 
Palomar 60-inch telescope performed photometry. RXTE PCA 
observed PKS 0537-286 on November 1 " 2006 in the 3 - 20 keV 
band. Details of the multiwavelength campaigns can be found in 
Tables [TJto 

2.1 .1 . UV to optical observations 

UVOT. - PKS 0537-286 was targeted in the U-b and with the 
UV-Optical Telescope UVOT (iRoming et alj|2005l) onboard the 
Swift satellite. Observations were performed in three epochs in 
2006 October and two epochs in 2008 February (Table [TJ. Here, 
we use the standard pipeline reduced imageproducts, co-added 
and exposure corrected within the XIMAGBj environment. 

Palomar. - B and V-b and observations with the robotic 
Palomar 60-inch telescope dCenko et alj|2006l) were earned out 
in eight nights between 2006 Oct 26.4 UT and 2006 Nov 2.9 
UT (Table [2]i. Data were reduced with standard IRAP0 routines 

1 See http://heasarc.gsfc.nasa.gov/docs/xanadu/ximage/ximage.html 

2 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association for Research 



dTodvlll993l) . Photometric calibration was performed relative to 
the USNO-B1 catalog which leads to a systematic contribution 
to the photometric uncertainties of ~ 0.5 mag. 

REM. - V, R, J, I, H observations were obt ained with the 
Rapid Eye Mount (REM) dChincarini et alj2003l) in October and 
November 2006 (Table [3). The REM data w ere reduced follow- 
ing standard procedures (Dolcini et al. 2 0051) . 

GROND. - The Gamma-Ray Burst Optical and Near 
Infrared Detector (GROND) is a 7-channels imager primarily 
designed for fast fo llow-up observations of GRB afterglows 
(Grei neret all 120081) . The use of dichroic beam splitters al- 
lows simultaneous imaging in seven bands, g'r'i'z' (similar to 
the Sloan system) and JHK (400 nm - 2310 nm). GROND is 
mounted on the 2.2 m ESO/MPI telescope on LaSilla/Chile since 
April 2007. The observations were done in 2 epochs in February 
2008 (Table©. 

2.2. X-ray observations 

INTEGRAL. - The INTEGRAL sat ellite observed the source 
for 1 Ms with its imager IBIS/ISGRI dLebrun et alj|2003l) which 
operates in the range of 17 - 1000 keV. The long exposure of 
1 Ms was carried out in non-consecutive pointings. These ob- 
servations were performed with a rectangular 5x5 dithering 
pattern and reported in T able [5l We used the sta ndard Off-line 
Science Analysis (OS A - ICourvoisier et al.ll2003l) software ver- 
sion 7.0 for the ISGRI analysis. Using the detected count rate 
we performed data cleaning. Most recent matrices available for 
standard software (isgri_arf_rsp_0025.fits) were used for spectral 
analysis. We used 296 science windows for a total amount of 98 1 

in Astronomy, Inc. under cooperative agreement with the National 
Science Foundation. 
3 See http://www.nofs.navy.mil/data/fchpix 
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MJD 




UTC 


Filter 


Mag 


Err 


54026.20225 


2006- 


10-18 04:51:14 


V 


18.25 


0.08 


54036.14096 


2006- 


10-28 03:22:58 


V 


18.83 


0.18 


54038.24992 


2006- 


10-30 05:59:53 


V 


18.50 


0.08 


54039.22512 


2006- 


10-31 05:24:10 


V 


18.46 


0.10 


54040.25737 


2006- 


11-01 06:10:36 


V 


18.39 


0.09 


54026.24084 


2006- 


10-18 05:46:48 


R 


17.67 


0.08 


54036.17956 


2006- 


10-28 04:18:33 


R 


17.93 


0.09 


54038.17376 


2006- 


10-30 04:10:12 


R 


18.36 


0.12 


54039.26192 


2006- 


10-31 06:17:09 


R 


18.55 


0.09 


54040.29142 


2006- 


11-01 06:59:38 


R 


18.48 


0.10 


54026.27942 


2006- 


10-18 06:42:21 


1 


17.15 


0.09 


54038.21236 


2006- 


10-30 05:05:47 


I 


17.74 


0.12 


54039.20966 


2006- 


10-31 05:01:54 


1 


17.77 


0.13 


54040.23736 


2006- 


11-01 05:41:47 


1 


17.77 


0.11 


54026.18431 


2006- 


10-18 04:25:24 


J 


17.05 


0.18 


54036.12300 


2006- 


10-28 02:57:07 


J 


17.68 


0.20 


54038.19724 


2006- 


10-30 04:44:01 


J 


17.66 


0.17 


54039.17244 


2006- 


10-31 04:08:18 


J 


17.65 


0.14 


54040.18487 


2006- 


11-01 04:26:12 


J 


17.73 


0.14 


54026.22602 


2006- 


10-18 05:25:28 


H 


16.43 


0.16 


54036.16471 


2006- 


10-28 03:57:10 


H 


16.77 


0.14 


54038.15901 


2006- 


10-30 03:48:58 


H 


17.41 


0.11 


54039.21766 


2006- 


10-31 05:13:25 


H 


17.10 


0.08 



Table 4. Log of GROND observations. 



Date (UTC) 




Filter 


Mag 


Err 


2008-02-14 02 


11:05 


g' 


19.25 


0.03 


2008-02-14 02 


11:05 


r' 


If 


3.89 


0.04 


2008-02-14 02 


11:05 


i' 


If 


3.86 


0.04 


2008-02-14 02 


11:05 


z' 


If 


3.88 


0.05 


2008-02-14 02 


11:05 


J 


If 


3.81 


0.06 


2008-02-14 02 


11:05 


H 


If 


3.76 


0.07 


2008-02-14 02 


11:05 


K 


If 


3.58 


0.06 


2008-02-14 05 


21:55 


g' 


19.25 


0.03 


2008-02-14 05 


21:55 


r' 


If 


3.89 


0.04 


2008-02-14 05 


21:55 


i' 


If 


3.85 


0.04 


2008-02-14 05 


21:55 


z' 


If 


3.82 


0.05 


2008-02-14 05 


21:55 


J 


If 


3.85 


0.06 


2008-02-14 05 


21:55 


H 


If 


3.75 


0.07 


2008-02-14 05 


21:55 


K 


If 


3.66 


0.06 



ksec of exposure time on the source. The source is detected with 
ISGRI in the 17-55 keV energy range with a significance of 5.6 
sigma. Data screening was performed according to the median 
count rate with respect to each science window and its distri- 
bution. Because of the low detection level a spectrum could not 
be extracted. The adopted dithering pattern resulted in a limited 
coverag e of the source b y JEM-X (Joint European Monitor for 
X rays - iLund et aLI2 003). For 77% of the time the source is out- 
side the FOV. For the remaining time the source is at angles < 4° 
from the axis. The total JEM-X useful exposure time is only 230 
ks, insufficient for a detection. 

Swift/BAT. - The Burst-Alert Telescope dBarthelmy et alJ 
120051) on board the Swift satellite is a coded-mask telescope op- 
erating in the 15 - 200 keV energy range. Thanks to its large 
field of view, BAT surveys up to 80% of the sky every day. We 
selected all observations of PKS 0537-286 comprised in the time 
span January 2006 - March 2008 . The data were pro cessed using 
HEAsoft 6.3 and following the recipes presented in lAjello et alJ 



Rev. 


Start date (UTC) 


End date (UTC) 


Exp (sec) 


0493 


2006-10-26 23:31:06 


2006-10-29 04:19:28 


180423 


0494 


2006-10-29 22:36:24 


2006-11-01 12:51:01 


202377 


0549 


2007-04-12 23:33:26 


2007-04-15 03:02:26 


170054 


0550 


2007-04-15 17:04:49 


2007-04-18 01:32:37 


188714 


0552 


2007-04-22 07:17:30 


2007-02-24 01:23:05 


136084 


0558 


2007-05-09 11:34:04 


2007-05-10 21:03:40 


115482 



Table 6. Log of X-ray observations. 



Instrument 


Start date (UTC) 


Exp (sec) 


XRT 


2006-10-27 05:55:53 


3000 


XRT 


2006-10-30 06:15:21 


3800 


XRT 


2006-10-31 06:21:26 


3400 


XRT 


2008-02-10 16:22:48 


6700 


XRT 


2008-02-12 16:24:38 


5200 


RXTE 


2006-11-01 00:06:40 


7200 


XMM 


2005-03-20 16:05:56 


80000 


BAT 


2006 - 2008 


2-years 



(l2008al) .The spectrum and the light curv e of PKS 0537-286 w ere 
extracted using the method presented in Aiel lo et ail ((2008b). 

RXTE. - The Rossi X-ray Timing Explorer (RXTE) ob- 
served PKS 0537-286 for 9.9 ksec, starting on 2006 November 1 . 
We reduced the data using the standard reduction script REX in- 
cluded in the HEAsoft6.0.4 package. The analysis was restricted 
to standard 2 binned data of layer 1 of the Proportional Counting 
Unit 2 (PCU2) of the Proportional Counter Array (PCA) from 3 
- 20 keV. We used PCARSP to produce the response matrix for 
the dataset. 

Swift/XRT. - The X-Ray Telescope (XRT) dBurrows et al.1 

12005b on board Swift observed the blazar in October 2006 and 
February 2008. The monitoring campaign consisted of five ob- 
servations. Data processing, screening and filtering were done 
using the FTOOL xrtpipeline included in the HEAsoft 6.3 distri- 
bution. 

3. XMM-Newton archive data 

We searched for observations of PKS 0537-286 in the XMM- 
Newton archive and found a 80 ks observation of the source 
(obs-id 0206350201) still unpublished. The source has been ob- 
served with the EPIC imagers. Data were linearized by applying 
the most recent calibrations with the software XMM-SAS ver- 
sion 9.0. We have extracted the light-curve of the background 
in the 10-15 keV energy range and excluded with a 3cr clipping 
technique a time interval of flaring particle background. Since 
the spectrum of many soft-proton flares may be extremely soft, 
we repeated the light-curve screening in the 0.3-10 keV energy 
band and excluded time interval with high background flux. As 
a result, most of the observing time has been lost because of par- 
ticle flares and the final sum of the good time interval yields an 
effective exposure of ~8 ks. 

By using the same software package we extracted the EPIC- 
PN spectrum of the source in a circle with a radius of 40 arcsec. 
The background has been evaluated on a similar area of the same 
PN-CCD, but far from source readout features. As a result the 
source spectrum contains ~3700 net counts in the 0.6-10 keV 
band which ensure one order of magnitude better statistics in 
comparison with the Swift data. 
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We further analyzed the XMM-Newton data of PKS 0537- 
286 observed during orbit 5 1 (but with la test response) in 199 9 
and reproduce the fit results as reported in lReeves et al.l d200ll) . 

4. Results 

4. 1 . Variability analysis 

Each single optical light curve obtained with REM, GROND 
and P60 does not show significant variations. Whenever, if at 
a given date, more than one REM, GROND or Palomar mea- 
surements are available, then the average flux is taken. Among 
the observations performed with the different optical instruments 
(this means among long time scales) the source showed variabil- 
ity as can be seen in Tables [2] [3] |4] Within each individual 
SwiftpCRT pointing just very marginal flux variation in the XRT 
spectrum was found (Table |7). Considering the large uncertain- 
ties among the five pointings, the flux amplitude is ~8%. The 
BAT light curve is consistent with being constant over the 2 two 
years of the survey. 

4.2. Swift and RXTE spectra 

The fit results of the single spectra of each single instrument are 
shown in Table [7] The spectra were fitted using XSPEC 12 and 
the latest available response matrices for calibration. The spectra 
are exceptionally hard. Swift/BAtT detected, in the 15 - 150 keV 
range, PKS 0537-286 with a flux of ~ 3 x 10~ n erg cirT 2 s" 1 
during the two y ears of the survey. Th e low Galactic N//-value 
(2.0 x 10 2l) cm- 2 lKalberlaetai] J2005h ) does not affect the BAT 
energy range, and we do not see evidence of further absorption 
or deviations from a single power law. 

Since the lower limit of the energy range is around 2 keV 
also the RXTE spectrum is not affected by absorption features. 
Therefore also this spectrum is well represented by a simple 
power law. 

The best fit results (x 2 red < 0.9) for each single SwiftfKRT 
spectrum is given by an absorbed power-law model, with ab- 
sorption parameter free to vary. Given the large fit uncertainties, 
reflecting the poor quality of these spectra, it is impossible to 
fully constrain the obtained parameters. 

Since we would like to assess the spectral characteristics 
with the best possible S/N, we have joined the average BAT, 
XRT and RXTE spectra and fitted the combined spectrum. The 
fit result are listed in Table [7] 

In the case of the co-added XRT, RXTE and BAT spectra 
we model the data with a power law with Galactic absorption 
controlled by the hydrogen column density. The large errors of 
the fitted parameters do not allow to constrain the absorption 
parameter. 

4.3. XMM spectra 

First we have checked the consistency of the single X-ray spec- 
tra among the 3 EPIC detectors within the sin g le ob servations 
by fitting power-law model as in iReeves et al.l (1200 ll) . Our re- 
analysis of the observation in 199 9 (orbit 51) confirms the re- 
sults reported bv lReeves et al.l d2001l) . where a simple power-law 
model fits the spectrum to a spectral index F = 1 .22 ± 0.08 with 
X 2 (dof) = 983(808). By fitting these data to an absorbed power 
law, the absorption parameter is consistent with the Galactic hy- 
drogen column density. Therefore no absorption in addition the 
the Galactic Nh is r equired, confirming once more the results in 
IReeves et alj d2001b . 



Table 7. X-ray fit results. 



Inst. 


obs epoch 


fit model 




r 


norm 


flux" 


(1) 


(2) 


(3) 


(4) 

n nol).26 


(5) 

1 


(6) 

T nfc.2 


(7) 

i a'I.X 



XRT 
XRT 
XRT 
XRT 
RXTE 
BAT 
joint 



2006-10-30 
2006-10-31 
2008-02-10 
2008-02-12 
2006-11-01 
2006 - 2008 



abs pi 
abs pi 
abs pi 
abs pi 

Pi 

Pi 
abs pi 



°0.00 

15 020 

005^ 

06" 18 
"•" u o.oo 



05° 08 

"•"-'0.02 



17 S5 

1 1:9 

1 5 1 9 
1 3 14 

• • J 1.2 



6 9 9:S 
3.6=' 

c nti 
J - '0.0 

14 7 ' 
4 2 4 6 



2. 

2.4* 

24 

2.3: 
6.3 
0.27^ 



"2.0 



i 

'3 



Explanation of columns: (l)=instruments; (2)=observation date; (3)=fit 
model; (4)=column density (units of 10 22 cm -2 ), (5)=spectral index; 



cm ); (7)=flux units are 



(6)=normalization (units of 10 4 ph keV~' 
10" 12 erg cnT 2 s for all instruments 

"fluxes are computed in the following energy ranges for the different 
instruments: XRT 0.6 - 6.0 keV , RXTE 3 - 20 keV, BAT 15 -150 keV, 
joint 0.6 - 150 keV. 




channel energy [keVJ 

Fig. 1. Data-to-model ratio for PKS 0537-286 as observed by 
XMM in 2005. The fitted model is a power law (including 
Galactic absorption) in the 0.3 - 8.0 keV energy range. The 
model is inadequate. The data are rebinned for demonstration 
to have at least a significant detection of 15 cr. The pn data are 
represented by squares, while circles and stars represent MOS-1 
and MOS-2, respectively. 



Contrary to the observation in 1999, the spectra derived from 
the observation in 2005 show an apparent absorption in excess 
to the Galactic value at ~1 cr. Figure Q] shows the data-to-model 
ratio of the best-fit power law modified by the Galactic absorp- 
tion (with ^ 2 (dof) = 1110(1013)). Besides falling to ~ 0.8 for 
the softest energies, the residuals exhibit emission in excess be- 
tween 1 -2 keV. This i s simila r to what is found by IYuan et all 
(120051) . IWorsley et all (l2004al) and iTavecchio et all (120071) for 
other high-reshift blazars. Therefore is the spectrum modeled 
by an absorbed power law with^ ,2 (dof) = 1015(1012). Figure 
|2] shows the confidence level contour plot of the absorption in 
excess to the Galactic value N# and the spectral index. In order 
to reproduce this clear spectral feature, we add another compo- 
nent to the base line model (wabs*powerlaw) and fit the photo- 
electric absorption component fixed to the Galactic value con- 
volved with the intrinsic absorption component and a power 
law (wabs(zwabs*powerlaw)). The resulting intrinsic absorption 
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Table 8. Fit results of X-ray XMM spectra. 



inst 


model 




r 


norm 


KT 


norm BB 




X Z (dof) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


PN-M1-M2 


ab*pl 


2.0 


1 24'-^ 


4 5 4 - 6 
50** 








1110(1013) 


PN-M1-M2 


ab*pl 


4 5 51 








1050(1012) 


PN-M1-M2 


ab(ab-*pl) 


2.0 


1-2211 


4 8 3 ' 






4 8 59 


1047(1012) 


PN-M1-M2 


ab(pl+bb) 


2.0 


l.OSjf 


t q4.0 


35° 40 






1057(1011) 



Explanation of columns: (l)=instruments; (2)=fit model; (3)=column density (units of 10 20 cm -2 ), (4)=spectral index; (5)=normalization (units 
of 10~ 4 ph keV -1 cm -2 ); (6)=photon temperature; (7)=blackbody normalization; (8)=absorption parameter at the source (units of 10 21 cm -2 ); 
(9)=statistical fit result. 




Fig. 2. Confidence level contour plot of the spectral index and the 
absorption parameter at 68%, 90% and 99% confidence levels of 
the spectrum extracted from the XMM observation of 2005. 



value is Ni 



4.8xl0 21 cm -2 . 



Tabel [8] shows the fit results of 
the co-added PN, MOS-1 and MOS-2 spectrum. 

A fit of the spectrum with a power law and black body, mod - 
ified by Galactic absorption as suggested in lCelotti et al.l d2007). 
gives an acceptable fit (x 2 red = 1.04). But there is only marginal 
improvement of the x 2 e( t (1.09) obtained fitting the power law 
with absorption fixed to the Galactic value. We further show 
that the same spectrum is modeled with x 1 , - 1 .04 by a power 
law and free absorption parameter. Figure |3]shows the unfolded 
spectrum fitted to the latter model in E 2 f(E) space. 

A systematic spectral flattening of high- redshift blazars at 
low X-ray energies wa s found pr eviously ( Fiore et al.l 1 1 9981: 
iReeves & Turner! 120001) . Similarly IWorslev et all d2004bl) de- 
tected a convex spectrum for the blazar PMN J0525-3343 at z 
= 4.4, and argued that it is most likely produced by a warm (ion- 
ized) absorber intrinsic to the source. We probed this hypothesis 
for PKS 0537-286 by fitting a power law and the contribution of 
ionized media intervening between source and the observer that 
modify the source flux in an energy-dependent wayQ. The warm 
absorber hypothesis results in an absorbing ionization state £ = 
150* 5 ° erg cm -2 s _1 . 

The XRT flux in the 0.6 - 6.0 keV band is systematically 
lower than the flux in the same energy band seen in XMM. 
Moreover, the spectral index of the XMM spectrum is slightly 
harder than the one of the XRT spectra. This suggests significant 




Energy (keV) 



4 The model is a convolution of three XSPEC models, namely pow- 
erlaw, absori and wabs fixed to the Galactic value. 



Fig. 3. Unfolded XMM spectrum between 0.3 and 8 keV in the 
vF v space. The spectrum is modeled with an power law and ab- 
sorption parameter free to vary. 



X-ray flux and spectral variability between 2005 and 2006-2008, 
so that we have not attempted any joint fit of the RXTE, Swift and 
XMM spectra. 

5. The spectral energy distribution 

We constructed two multiwavelength spectra using the data of 
the first and fourth epochs (corresponding to the XRT pointings 
on October 27 f/l 2006 and on February 10 ,h 2008) of our cam- 
paign, and adding to these non-simultaneous multiwavelength 
data, including an EGRET detection in three energy bins in the 
GeV range. The optical and UV data w ere corrected for Ga lactic 
absorption using E(B-V) = 0.025 mag dSchlegel et al.l 19981) . and 
converted into fluxes. The optical flux at wavelengths shorter 
than ~5000 A is significantly suppre ssed by the Ly a interven- 
ing absorption ( We vmann et al.lfl98 lh . Therefore is the sharp in- 
crease from UV towards optical band not surprising, as the Ly a 
forest is shifted along lower frequencies due to the high-redshift. 
PKS 0537-286 belongs to the FSRQ subclass of blazars. The 
typical two hump structure is evident in the two SEDs. The IC- 
peak can be constrained by the BAT data points and EGRET data 
to be located at a few GeV. The synchrotron peak is somewhat 
difficult to determine. 

We model the data according to iGhisellini & Tavecchiol 
(2009), assuming a one-zone, homogeneous synchrotron and 
inverse-Compton model. The ejecta move within the jet seen un- 
der an angle and bulk Lorentz factor T. The homogeneous mag- 



6 E. Bottacini et al.: PKS 0537-286, a multifrequency campaign 



Table 9. List of parameters used to construct the theoretical SED. 



^diss 


M 


^BLR 


F 


u 


B 


r 


By 


7b 


ymax 


Si s 2 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) (12) 


4.2 (2) 


7 


1.3xl0 J 


0.3 


168 (0.16) 


0.1 


18 


3 


40 


10 4 


-1 2.6 


3.2(1.5) 


7 


916 


0.2 


84 (0.08) 


0.09 


18 


3 


25 


10 4 


-1 2.6 



Explanation of columns: (1) = dissipation distance in units of 10 18 cm and (in parenthesis) in units of 10 3 R$\ the (spherical) emitting size is ten 
times smaller; (2) = black hole mass in 10 9 solar masses; (3) = size of the BLR in units of 10 15 cm; (4) = power injected in the blob calculated 
in the comoving frame, in units of 10 45 erg s ; (5) = accretion-disk luminosity in units of 10 45 erg s -1 and (in parenthesis) in units of LEdd; (6) = 
magnetic field in Gauss; (7) = bulk Lorentz factor at ^dj ss ; (8) = viewing angle in degrees; (9) and (10) = break and maximum random Lorentz 
factors of the injected electrons; (11) and (12) = slopes of the injected electron distribution (<2(y)) below and above yb- 



netic field is tangled throughout the emitting region. This is lo- 
cated at a distance R^ ss from the black hole, has a cross sectional 
radius R = O.l^diss and a (comoving) width AR' = R. The par- 
ticle distribution is derived through a continuity equation, and 
assume that particles are injected for a finite time, of duration 
tj„j = AR/c. Radiative cooling and electron-positron pair pro- 
duction and emission (not relevant in this case) are accounted 
for. Although the considered continuity equation has no steady 
state solution, we evaluate it at the time ti„j, when it reaches the 
maximum value, corresponding to the maximum emitted lumi- 
nosity. The broad line region (BLR) clouds reprocess and re- 
isotropize 10% of the disk luminosity. Since we assume that the 

1 /2 

radius of the BLR is oc L ' , the radiation energy density of the 
broad lines (within the BLR) as seen in the comoving frame de- 
pends only on T 2 : t/ BLR oc T 2 Lblr/^blr =constxF 2 . This com- 
ponent, as seen in the comoving frame, is approximated with a 
black body spectrum (see Tavecchio & Ghisellini 2008). We as- 
sume that the injected particle distribution is a smoothly broken 
power law, with slopes s\ and S2 below and above the electron 
Lorentz factor y^. In the case of our source, radiative cooling is 
efficient, implying that electrons of all energies have time to cool 
in ti„j. As a consequence, the particular value of si is ininfluent, 
since, below y\,, the emitting particle distribution is oc y~ 2 , while 
it becomes oc above y\,. We modeled the first and the 

fourth observing epochs which show the largest variation from 
near-IR to UV band. Using the detection in three energy bins 
by EGRET, the IC peak was assumed to be between these data 
points and the Swift/BAT spectrum. The results are shown in 
Figure |4](solid orange and blue lines). The assumed parameters 
are reported in Table [9] 

Table [10] lists the power carried by the jet in the form of 
radiation (P r ), magnetic field (Pb), electrons (P e ) and protons 
(P p , assuming one proton per emitting electron). The powers are 
calculated as 

Pi = nR 2 r 2 /3cUi (1) 

where {/, is the energy density of the i th component. 

The optical black data points are not well represented by the 
model. Indeed, assuming the highest data point to be the peak of 
the accretion-disk component, then the model requires a factor 
~3 more in temperature. Being the relation between accretion- 
disk luminosity (Lj), accretion rate (M) and temperature (T): 

oc M oc T 4 , (2) 

that, in turn would require a very high variation of the accretion 
rate (factor ~80), at odds with the model. Therefore we decided 
to assume the peak of the accretion-disk component around the 
GROND z' data point, where a turn over in the data sequence 



can be seen and that can well reproduce the IC component. The 
jet power and the observed non-thermal luminosity increase as 
the disk luminosity increases, despite the fact that the amount of 
inverse Compton cooling is constant (i.e. U' BLR is constant). The 
variations of the non-thermal bolometric luminosity are due to 
the different power injected in the form of relativistic electrons, 
and the different Compton to synchrotron luminosity ratio is due 
to the different assumed magnetic field. 



Table 10. List of parameters used to construct the theoretical 
SED. 



logP r 


logP B 


logPe 


logPp 


(1) 


(2) 


(3) 


(4) 


46.5 


45.3 


46.6 


48.0 


46.2 


45.0 


46.4 


48.0 



Explanation of columns: (1), (2), (3) and (4): jet power in the form of ra- 
diation, magnetic field, bulk motion of electrons and protons (assuming 
one proton per emitting electron). These powers are derived quantities, 
not input parameters. For all cases the X-ray corona luminosity L x = 
0.3Ld- Its spectral shape is assumed to be oc v -1 exp(-/;v/150 keV). 



6. Discussion 

It is unlikely that the apparent excess column density for PKS 
0537-286 is of Galactic origin. It would require absorption by 
cold material with a column of N# ~ 2.5xl0 20 cirT 2 in ad- 
dition to the Galactic level toward the direction of the source. 
The Galactic hydrogen colu mn density derived wit h in different 
energy ban ds (IR, radio) by ISchlegel et alj d 1998b . lElvis et al.1 
d 19891) and iKalberla et all d2005l) enforces the accuracy of the 
measured Galactic N#. Also the intergalactic origin is highly im- 
probable because the line of sight would have to intercept several 
systems with very high column density. The chance probability 
to intercept several very high col umn density Ly a absorptio n 
systems is very low (p ~ 2%) (O'Flaherty & Jakobsen 1997). 
Furthermore, we exclude that the damped Ly a absorption sys- 
tems (DLA) contribute significantly to the spectral absorption 
feature. In fact, the DLA at z = 2.974, associated with an inter- 
vening galaxy, has HI column density of logNni = 20.3 + 0.1 
dAkerman et al.l 120051) and metallicity below half solar. This 
would give an absorption at soft X-ray energies caused by a 
column density of N# < 10 20 cm -2 . The intrinsic absorber hy- 
pothesis is much debated. For this absorber we derived a high 
column density value of ~ 4.8xl0 21 crrT 2 . Under the warm 
absorber hypotesis we cannot constrain the ionization state of 
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Fig. 4. Observed spectral energy distribution of PKS 0537-286 on 27 October 2006 (black filled circles) and 10 February 2008 (red 
filled circles). The filled green triangles represent the Swift/BAZT spectrum extracted from the 2 years survey. The data are corrected 
for Galactic extinction. Gray points are archival data. The dotted line shows the single components to the SED from the IR torus, 
the dis k and the X-ray corona . The p urple solid line represents the 5 <x sensitivity of Fermi after one year mission time. The jet 
model Ghisell ini & Tavecchiol (120091) for the energy distribution for 27 October 2006 (orange solid curve) and 10 February 2008 
(blue solid curve) are overplotted. 



the possible absorbing material better than £ = 150q q erg cm 
s _1 . Furthermore, the best evidence for warm absorbers, like the 
O absorption edge and the Fe M-shell, are down shifted outside 
our observed range (0.3 - 10 keV) being in the rest frame's en- 
ergy band within 0.5 - 0.9 keV. The possibilit y for the warm 
absor ber to give rise to absorption variability dReeves & Turner! 
2000) is not confirmed in our case. 

The excess absorption, also referred to as spe ctral flatten- 
ing, has been reported several times for this so urce (ICappi et al.l 
119971: iFiore e"taLlll998t ISambruna et alj|2007l) . Spectral flatten- 
ing is in agreement with a general behavior seen in a number 



of high-redshift radi o-loud quasars studies (IFiore et al 1 fl998l 
iReeves & Turn er 2000), where the required column density is 
up to \0 22 cmr 2 , whereas low-redshift blazars show unabsorbed 
spectra. If the flattening is not a mere selection effect due to the 
high-redshift, then this could indicate an ev olutionary scenario. 
In fa ct, the positive evolution of FSRQs dDunlop & Peacockl 
119901) foresees a bulk of FSRQs at high-redshift. PKS 0537- 
286 shows time variability of the apparent excess column den- 
sity. Figure [5] shows the time variability of the column den- 
sity over 7 ye a rs as measured by several instruments. Similarly 
Scha rtel et al.l (1 19971) find variable absorption of X-rays for a 
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- 2 Lj i i i i i Ll: 

1990 1995 2000 2005 

UTC (years) 

Fig. 5. Time variability of Nh measured by several instruments. 
From left to right the instrume nts are: ROSAT (ICappi et alj 
Il997l) . two ASCA upper limits (ISiebertet al.ll 19961) . XMM (this 
work), Chandra 3 cr upper limit(this work), XMM (this work). 
The dashed horizontal line delim its the Galactic hydrogen col- 
umn density dKalberla et al.ll2005l) . 

high-redshift blazar on time scales of < 1.5 years. At present 
the origin of this absorption is unknown. 

Inspection of the SEDs allows us to constrain the inverse 
Compton peak rather well, thanks to the exceptionally hard X- 
ray spectrum and the EGRET data. Moreover, o ur source is not 
included in the Fermi bright source list (lAbdo et al.ll2009l) . in- 
dicating that the upper limit is consistent. The model does not 
account very well for the shorter wavelength optical data, be- 
cause of the Ly a drop. An optical flux excess in the 2 nd epoch is 
partially due to strong emission lines contaminating the GROND 
filters. The chosen parameters for the I st observ ation epoch de- 
viate f rom those obtained for the same source bv | Sambruna et aTl 
d2007l) and modeled using lGhisellini et al.l (120021) . This is due to 
our accurately sampled data, that allow us to better constrain the 
theoretical SED. 

7. Conclusions 

We have conducted a multifrequency campaign on PKS 0537- 
286 on six different epochs. The source was monitored nearly 
simultaneously from optical to gamma-rays. In the optical - near- 
IR band the most evident flux variations were observed between 
the observations on October 27' /! 2006 and on February 10 th 
2008. The compiled SED for the latter two observation epochs 
allowed to derive the physical parameters of the source. 

Also a new and detailed X-ray spectroscopic study of the 
XMM-Newton observation in 2005 of the source has been pre- 
sented. The spectral flattening for this sour c e at so f t X-ray en- 
ergies , repo rted previously by Cappi et al.l d 1997b . iFiore et all 
d 1998b and ISambruna etall (120071) . is confirmed by the high 
signal-to-noise XMM spectra. This flux decrease at soft X-ray 
energies cannot be attributed to excess absorption in the Milky 
Way or to intervening material in the intergalactic medium along 
the line of sight to the source. Data suggest that the absorber 
can be intrinsic to the source. By using our data and search- 
ing in literature we have constructed a time curve of the Nh 
for the source. We confirm the variability of the apparent excess 
absorption on time scale of years. A periodicity for this spec- 
tral feature cannot be inferred. Our result highlights once more 



the evidence for the existence of a population of superluminous 
blazars at high-redshift that show a soft X-ray cut-off. If the bulk 
at high-redshift of these type of sources is not only due to se- 
lection effect, then their existence i s important in the context of 
cosmological evolution of SMBHs (Bottc her & Derm er 2002). 
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